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Abstract 

We study the finite term of the holographic entanglement entropy for the charged black hole 
in AdSd+2 and other examples of black holes when the spatial region in the boundary theory 
is given by one or two parallel strips. For one large strip it scales like the width of the strip. 
The divergent term of its expansion as the turning point of the minimal surface approaches the 
horizon is determined by the near horizon geometry. Examples involving a Lifshitz scaling are 
also considered. For two equal strips in the boundary we study the transition of the mutual 
information given by the holographic prescription. In the case of the charged black hole, when 
the width of the strips becomes large this transition provides a characteristic finite distance 
depending on the temperature. 



Introduction 



Entanglement entropy is an important quantity which has been studied in many models of 
condensed matter systems, quantum information and quantum gravity. It measures the quantum 
correlations in a bipartite decomposition of a quantum system. 

Let us consider a system whose total Hilbert space can be written as a direct product H = 
Ha ® Hb- Denoting by p the density matrix characterizing the state of the system, the reduced 
density matrix associated to A is obtained by tracing p over the degrees of freedom of B, 
i.e. pa = Tr^/j. Then, the entanglement entropy is defined as the corresponding Von Neumann 
entropy, namely Sa = — Tta(pa log Pa)- When the system is in a pure state, we have p = | V I / )( V I / | 
and Sa = Sb- A very interesting situation occurs when A and B correspond to a spatial 
bipartition of the system. In this case the entanglement entropy is called geometric entropy. 
Here we consider this quantity and we will always refer to it as the entanglement entropy. Its 
interesting feature is that is satisfies the so called area law: the leading term in the expansion for 
small UV cutoff a is proportional to the area of the boundary separating A and B. In d spatial 
dimensions we have Sa oc Aiea,(dA) / a d ~ l + . . . , where the dots represent higher order terms in 
a pp. This area law is violated in two dimensional conformal field theories, where a logarithmic 
behavior has been found for one interval. In particular Sa = ( c /3) log(^/a) where i is the length 
of the interval and c is the central charge of the theory. The method employed to get the analytic 
result for Sa is the replica trick, which means first to compute Trp^ for integer n and then to 
perform an analytic continuation to real values of n in order to take Sa = — 9 n Tr/9^|| ri=1 [2JE1I1] 
(see [5] for a recent review). 

For quantum field theories with a holographic dual, the problem of computing the entanglement 
entropy through a bulk description has been addressed in [61 [7] . The holographic prescription 
to obtain Sa associated to a region A in the d+1 dimensional boundary theory is the following. 
On a fixed time slice (see [8] for a generalization to time dependent backgrounds), among all 
the d dimensional surfaces extended in the bulk whose boundary coincides with the boundary 
of A, we have to consider the one having minimal area. Denoting this minimal surface by 7,4, 
the holographic entanglement entropy is given by S A = Area( 7 A)/(4G^ +2) ), where G$ +2) is the 
Newton constant of the d + 2 dimensional theory in the bulk. Besides recovering the area law, 
this prescription passed many tests (e.g. the strong subadditivity inequalities) and it has been 
deeply studied (see the recent review [9] and the references therein) ; thus it is considered a key 
tool to understand the essential features of the entanglement entropy for quantum field theories 
with a holographic dual. 

The entanglement entropy is not an extensive quantity, as can be easily understood e.g. by 
the fact that 5,4 = Sb for a pure state (this equality is violated at finite temperature). In the 
holographic computation of the entanglement entropy extensivity is recovered if one considers 
the finite term of the minimal area (sometimes called renormalized entanglement entropy), i.e. 
the one obtained by subtracting the UV divergent term giving the area law, in a black hole 
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background [TJ ITUl [TTj I12j . This behavior is due to the fact that, as the size of the region A in 
the boundary tends to infinity, a large part of the minimal surface gets very close to the horizon 
and this part goes like the volume of A in the large size limit. Thus the near horizon geometry is 
responsible of the leading divergence of the finite term of the minimal area as the turning point 
of the minimal surface approaches the horizon. 

A second important aspect of the entanglement entropy we are interested in concerns the case of a 
spatial region A in the boundary made by two disjoint regions, i.e. A = A1UA2 with A\C\A2 = 0. 
In this case the natural quantity to consider is the mutual information M4 = Sai + Sa 2 — Sa 1 uA 2 
because it is UV finite. For some spin chains and two dimensional conformal field theories 
interesting results have been obtained jT3[ HH US HSJ H7J HHJ HH [20] . The models considered in 
these papers have small central charges (order of the unity) . 

In the context of the holographic correspondence, the case of disjoint regions has been addressed 
in [HI [71 [2"T] [2"2"[ [2U] I24j . An interesting feature of the holographic entanglement entropy is the 
transition of the mutual information from zero value to a positive value (the mutual information 
cannot be negative, as a consequence a strong subadditivity inequality) |21[ I24j. This transition 
should be a large c effect, which is the regime where the holographic prescription works, since 
there no signal of it e.g. for the compactified boson |14j. 115]. which has c = 1. 
In this paper we consider the holographic entanglement entropy for one or two strips in the 
boundary theory in presence of various types of black holes with non compact horizon in the 
bulk. For one strip, we focus on the divergence of the finite term when the strip becomes large 
and therefore the turning point of the minimal surface approaches the horizon. The degree of 
this divergence depends on the near horizon geometry, but the finite term scales like the width 
(and thus like the volume) of the strip for all the black holes we consider. This scaling is broken 
for the Lifshitz type backgrounds whose dynamical exponent occurs in the spatial part of the 
metric. 

For two parallel strips of equal width we mainly consider the transition of the mutual information 
in terms of the geometrical parameters, namely the width of the strips and the distance between 
them. For the charged black hole in four dimensions with fixed charge, we find that the transition 
of the mutual information leads to a characteristic finite distance between the strips as they 
become large. This distance depends on the temperature and it could be interpreted as a signal 
of the occurrence of a finite correlation length in the boundary theory. 

The paper is organized as follows. In the section [T] we review the holographic prescription 
for the entanglement entropy, specializing the analysis to ansatz that contain the black hole 
metrics we consider in the rest of the paper. In the section [2] we study the finite term of the 
holographic entanglement entropy for the charged black hole, the warped black hole of [25] and 
the perturbed Lifshitz background considered in [26 as a solution of the Abelian Higgs model 
[27] , In the section [3] we study the Lifshitz black hole of [28] computing the analytic expression 
of the holographic entanglement entropy to all orders in the UV cutoff. This allows us to extract 
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the finite term and to test the method employed for the other black holes. In the section [4] we 
consider two equal and parallel strips in the boundary and study the transition of the mutual 
information for AdSd+2 and for the charged black hole. 



1 Holographic entanglement entropy for black holes 

In this section we review the holographic prescription to compute the entanglement entropy 
O [7], defining the integrals we need in order to study the black holes that we will consider in 
the remaining sections. In the appendix |A| we review the results for AdSd+2, that will be also 
employed in the section |4j 

Given a quantum field theory living on the boundary Mf x R d of an asymptotically AdSd+2 space, 
we take a d dimensional region A strictly included in the constant time slice of the boundary. 
Let us take a <i-dimensional surface 7 embedded in the constant time slice of AdSd+2 defined by 
z = z(x), being z is the holographic coordinate and x 6 M. d a vector of the constant time section 
of the boundary. The area of 7 reads 



Area(7) = J dx\ . . . dxd \J det (hi j) (!•!) 



where hy is the induced metric on the surface dsf nd = hijdx l dx^ . Among all these surfaces, we 
restrict our attention to the ones whose boundary coincides with the boundary of the region A. 
Within this smaller subset of surfaces, we denote by 7^ the one having minimal area. 
The proposal of Ryu and Takayanagi [6J, [7] is that we can holographically compute the entan- 
glement entropy in the boundary theory through a computation in the bulk. In particular 

Area(7 A ) . , 

Sa = i5r>" <L2) 

where G^ +2 ^ is the Newton constant in d + 2 spacetime dimensions. 

Depending on the shape of dA, one decides if it is more convenient to work either in cartesian 
(dx 2 = ^f = i(dx 1 ) 2 ) or polar coordinates (dx 2 = dp 2 + p 2 dO,,_ 1 , being dVL 2 ,^ the metric of the 
d — 1 dimensional unit sphere) of M. d . Since we will mostly consider A to be a finite strip or a 
disjoint union of two of them, we will adopt the cartesian coordinates for W 1 (for an example 
where the polar coordinates system is employed, see the appendix [Cj which contains a discussion 
on the circular case in the black hole background considered below). 

For many known black holes which are asymptotically AdSd+2, the metric on the fixed time 
slice is given by 

d4^ds 2 l=R 2 ( d 4 + J^) (1.3) 



lt=const " ^ z 2 z 2 f(z) / 

where R is the radius of AdSd+2 realized closed to the boundary. In this system of coordinates 
the boundary is the z = slice and the horizon is characterized by the smallest zero of the 
emblacking function f(z). 
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Let us consider the region A in the boundary given by a strip with length L along one direction, 
that we call x, and L± along the other orthogonal ones. Choosing the origin in the center 
of this strip, the symmetry of the problem allows us to restrict to surfaces described by the 
even function z = z(x). Then, the area functional that we have to minimize to compute the 
holographic entanglement entropy reads 



Area(7A) 



2R d L d r l 




(1.4) 



Considering as a Lagrangian density £ s trip[z(aO] the integrand in (1.4), one notices that it does 
not depend explicitly on x. This is the main simplification that makes the case of a rectangular 
region A easier to solve than the case of a circular region. Indeed, the independence of £ s trip on 
x leads to the conserved quantity % s tiip = Pz%' — £ str ip, where p z = dC s t r i P /dz' . In particular, 



one gets % s tri P = z d [l + ( z ') 2 /f(z)] 1 ^ 2 - By introducing z 



2d 
max 



l/7ig tlip , the constancy of 



H 



strip 



reads 



z 



2d 



y 2d 



Z" 



where we have used that z' < 0. This equation tells us that z, 

z' = when z = z mSiX . Notice that, from (1.5), also at the horizon zq we could have z' 



(1.5) 

is the turning point, namely 





because f(zo) = 0, but we never reach it because zq > z mSLX ^ z(x) ^ 0. The equation (1.5) 
provides the profile of the minimal surface we are looking for and, by construction, it satisfies 
z{L/2) = and z(Q) 
from £ st rip 



Zmax- As a check of (|1.5|), one can write the equation of motion coming 

2 ft*) 



z"z 



d 



1 + 







(1.6) 



and verify that the same equation can be found by deriving the conservation law T~L 



strip 



const 



w.r.t. x. Then, separating the variables in (1.5), we find that the inverse function x(z) is 



x{z) 



dx 



w 



Vf( w ) V 41 - w 



,,2r/ 



: dlJJ . 



(1.7) 



Imposing in (1.7) the relation x(0) = L/2, one gets that 



L 

2 



w 



V7H V^max - W 



,,2d 



: dw 



(1. 



L(z max ) and the correspondence between z max and L. 



which provides L 

As for the area of the minimal surface defined by (1.5), we can employ its definition to change 



integration variable in (1.4), which therefore becomes 



Area(7A) 



2R d L d r x 



Z d y/f(z)y/z^~^Z 



>2d 



dz . 



(1.9) 
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It is important to remark that, since f(z) — > 1 as z — > 0, the integral in (1.9) diverges at z = 0. 



This leads us to put a UV cutoff z ^ a in the integration domain of (1.9). Thus, the integral we 
have to compute reads 

lax ~d 



Area(7 j4 ) = 2R a L 



d t d—1 



w d y / f(w)^z^ ax -w 



■lit 



dw 



R d L^ l A d {z m ^a) 



(1.10) 



In order to isolate the divergence of (1.10) as a — > 0, we write the integral as follows 



A d (z n 



ur 



dw + 



2,, 



wd Vy7HV4 d ax-^ 



2d 



1 



(d - 1) a d - 



(1.11) 



(1.12) 



The second integral in (1.11) is finite as a — > because we have either f(w) = 1 + 0(w d+1 ) or 
f(w) = 1 + 0(w d ) for w (see (2.3) and (3.1) respectively). In (1.12) we have introduced 



the finite term in the UV cutoff expansion 



wd \y/7MV*%*-™ 



2d 



1 dw 



(d-l)z 



d-l 

max 



;i.i3) 



In this paper we will be mainly interested O(l) term in the a expansion, which is ^(z max , 0). In 
the figure [T] this term is shown for the charged black hole in AdS^ at zero temperature. As the 
turning point z max approaches the horizon zq, it develops a divergence we are going to study. 
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Figure 1: Charged black hole, extremal case and zq = 1. Plot of the finite term ^2(z max , 0) as 
function of z max (red curve). Close to the boundary (i.e. when z max —¥ 0) it coincides with the 
curve corresponding to AdS^ (blue curve), which can be read from ( A.8| ). 
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Equivalently, we can isolate the divergence for small a in the integral of (1.10) as follows 

rZm&x o f2 max n / d \ 

A d (z max , a) = I - ,== (fa, + / - ,== ( , ^ - -lU (1.14) 



'Vfi 



w d y/f(w) I \Amax - w 



'Vfl 



in 



dw + I d (a, Zr, 



(1.15) 



The finite term in the expansion for small a is now given by I d (0, z ma , x ) plus a contribution from 
the first integral in ( 1.15| ). The distinction between ( |1.11[ ) and (1.14) is obviously meaningless 
for AdSd+2, where f(z) = l identically. The splitting ( |1.14[ ) has been used in the appendix [P| to 
get some insights about the expansion of the finite term of the minimal area as a power series 
as zq — z max and the possibility to approximate it through the near horizon geometry. 



1.1 A more general ansatz 

In this section we consider a more complicated expression for the metric in order to include 
other kind of black holes in our discussion. Let us take a D + 1 dimensional spacetime and the 
following ansatz for the metric on the constant time slice 

+ B(r) 2 dx 2 ) e-TxM + R 2 e^^? (r)d9W 
A{r) z J J 



dsn 



;i.i6) 



where dx 2 gives the metric of M. d and 7^ (r) is the metric of a D — 1 — d dimensional compact 
manifold M c . The boundary is at large r and we assume the occurrence of a horizon at r = r^. 
Let us take a strip specified by the function r = r(xd). Then, the metric induced on it reads 



ds 2 

as md 



B{r)\dx{ + ...dx 2 d _ 1 ) + 



Air) 



+ B(r) 



dx\ I ( 



(1.17) 

To compute the area, we have to integrate \J det(Gj n d) over the strip A. In such determinant 
the dependence on x( r ) simplifies and therefore it does not occur anymore. If det(7^ c ^) does not 
depend on r, then the area of the surface is given by 



Area(7A) = [R 




„A2 



dx d 



10 V A{r) 2 B{r) 2 

where L± is the width of the strip along the directions x±, . . . , Xd-i and Vol(At c )- 



(1.18) 



As done above, we take as Lagrangian density >C s t r ip[^(^d)] the integrand of (1.18) and compute 
the momentum p r = d£ s t r i p /dr' . Then, being £ s t r i P independent of x d , we can employ the 
conserved quantity % s trip = p r r' — jC s trip • Since at the minimum value of r we have r'(0) = 0, 



we set %s trip = B(r m i n ) 2d . This allows us to write (1.18) as follows 



Area( 7 A) = [R D ^ d \o\{M c )]2L 



d-l 



B{r) 



2d-l 



A(r)^B(r) 2d -B(r min ) 2d 



dr . 



(1.19) 
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It is also important to express L = L(r m i n ) and it reads 

-S(^min) 



L 



lin 4(r)5(r)VB(r) M -B(r mi 
We require to have AdSd+2 at large r, which means to impose 



\2d 



: dr 



(1.20) 



A(r) 



R 2 



O(l) 



i? 2 



+00 . 



(1.21) 



Because of this asymptotic behavior, the integral in (1.19) is divergent. Thus, one introduces 
the cut off a at large r, obtaining for the regularized area 

Area( 7A ) ord _x f a B(r) 2d - 1 



RD-i-dy ol ( Mc 



2L° 



pet 

1 1 



dr 



2L° 



(d - l)R d ~ 2 



a 



d-l 



„d-l 



in A(r)^B(r) 2d -B(r min ) 2d 



+ 2L 



d-2 



(1.22) 

dr 
(1.23) 



where the integral in (1.23) is finite when a — > 00, once the asymptotic behavior (1.21) has been 
assumed. At this point, the finite term of area integral we are interested in is given by the sum 
of the integral and of the term proportional to rf^^ in (1.23). We remark that (1.8) and (1.10) 



are special cases of (|1.20|) and (jl . 22|) respectively. Indeed they are recovered by choosing 

vW) 



A{r) 



r 
R 



B(r) 



r 
R 



(1.24) 



and adopting the variable z = R 2 jr. The formula for the holographic entanglement entropy 
then gives 



S A 



Area(7^) 



2L° 



AG 



(D+i) 



N 



AG 



{d+2) 
N 



B(r) 



2d-l 



lin A(r)^/B(r) 2d -B(r min ) 2d 



dr 



(1.25) 



where we have used that G^ +1 ^ = G^*' [R u ~ 1 ~ d Vol(A^ c )]- Notice that the compact part 
enters through Kaluza-Klein reduction in the Newton's constant also in this case where a warping 
factor occurs between the compact and the non compact part [11| IT2] . 



i{d+2),j-,E>-l-d~ 



2 Expansion of the finite term near the horizon 

In this section we study the finite term of the holographic entanglement entropy introduced in 
the previous section. In particular, we consider the leading term of its expansion as the turning 
point of the minimal surface approaches the horizon, which means that the width L of the strip 
in the boundary becomes large. As examples, we analyze the charged black hole in AdSd+2 



(section 2.1), the warped black hole of [25] (section 2.2) and the perturbation of the Lifshitz 



background found in [26] within the context of the Abelian Higgs model of |27j (section 2.3). 



The finite term in the expansion for small UV cutoff a is given by Ad(z max , 0), defined in ( jl. 13 ) . 
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In order to consider its expansion as the turning point z max of the minimal surface gets close to 
the horizon, we take 

^max = ZQ- eCmax £ ~> finite (max (2.1) 



and change the integration variable in (1.13) according to this expansion, i.e. we set w = zq — e£, 
where < Cmax < C- Then, the finite term can be written as follows 

P Cmax 

(d-l)z 



- '• '' ' • •" i::ax 



where £? C [/c m ; n ,oo) C Q is some discrete set of increasing rational numbers, which are not 
necessarily positive (fe mm < 0). For instance, in the case of the charged black hole with d = 2 
we have k £ { — 1/2, 1/2, 1,3/2, . . . }. In order to write Ad(z max ,0) as an expansion in terms of 
powers of zq — z max , we have to compute the definite integrals occurring at each k and then 
expand each of them for small s. Then this expansion can be written in powers of zq — z max by 



using the definition e = (zq — z max )/( max from (2.1). 



In all examples we have considered we find that this method provides only the divergent term 



as Zmax - >■ z o- This is due to the fact that all the integrals occurring in (2.2) give a contribution 
to the finite term of the expansion. 

The same procedure just described to expand the integral Ad(z ma , x ,0) can be applied to the 



integral (1.8) as well, obtaining L as an expansion in powers of zq — z max - It is then useful to 
compare the divergences of these two quantities as z max — > z$ in order to see how the finite term 
of the entanglement entropy scales with the width of the strip, and therefore with the volume. 

2.1 Charged black hole 

In this section we apply the method just described to the charged black hole in AdSd+2 in 
its three different regimes of neutrality, extremality and non extremality. The metric and its 
properties are reviewed in the appendix [Bj 
The metric of the charged black hole in AdSd+2 reads 

ds 2 -fdt 2 + dx 2 dz 2 „ ~ / z \^ %r / z \<*+i , . 

where M is the mass and Q is the charge of the black hole. The radial direction is parameterized 
by z and the boundary is at z = 0. The position zq of the horizon is given by the smallest zero 



of the emblacking function f(z). Since the metric (2.3) falls into the class of metrics described 



by (1.3), we can employ the formulas discussed in the section [TJ 



Schwarzschild black hole. As a first example, we consider the Schwarzschild black hole, 



which is given by (2.3) with Q = 0. By performing the expansion described above, we find 



L = — log(z " ^rnax) + 0(1) (2.4) 

Vd{d + 1) 
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and 



•^d(-2max) 0) 



V2_ 



^ log(z - Zmax) + 0(1) 



L 



+ 0(1) 



(2.5) 



where we recall that the horizon zq is related to the temperature as T = (d + 1)/(4tt2o)- The 
case d = 3 was considered in [7j . 

• Extremal charged black hole. When Q / and T = this analysis leads to 



L 



d^/d+l yfz^ 



and 



•^d(^mam 0) 



V2 



TT Z. 



3/2-d 



+ o(i) 



+ 0(1) 



L 



+ 0(1). 



(2.6) 



(2.7) 



+ 1 V z - 2max 

In the figure [2] (see |llj ) we show ^4rf(z max ,0) in terms of L for the extremal case. When L is 
small we are close to the boundary and the curve reproduces the one of AdS±, as expected. By 



comparing the two plots in the figure, one can check the dependence on zq in (2.7). 
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Figure 2: Extremal charged black hole in AdS^ with zq = 1 (left) and zq = 1.5 (right). Plot of 
the finite term ^(^maxi 0) as a function of L (red line). For small L it recovers the corresponding 



quantity for AdS± (blue curve) obtained from (A.9). The black line provides the large L behavior 



given by (2.7). 



Non extremal charged black hole. The same method applied for T > leads to 



L 



y/2irdT 



log(^ 



and 



*A-d(z max , 0) 



log(2 



+ 0(1) 



+ o(i) 



L 



+ 0(1). 



(2.8) 



(2.9) 



Comparing these three regimes of the same black hole, one learns that the finite term of the 
holographic entanglement entropy diverges like the width L (and therefore like the volume) of the 
strip in the boundary. The distinguished feature is the kind of divergence of Ad(z max , 0) and L as 
Zmax — > zq. This is determined by the near horizon geometry which is given by f(z) = 0(zq — z) 
for the Schwarzschild and the non extremal case and by f(z) = 0((zq — z) 2 ) for the extremal 
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case (see the appendix [B|) . As a check, one can perform the expansion of the finite term just 
described substituting to the emblacking function f(z) its near horizon approximation f n h(z) 
and verify that the same divergence shown above are obtained. 

We remark that for all the black holes we are considering the horizon is non compact; therefore 
the wrapping of the minimal surface around the horizon in the large L limit described e.g. in 



[3 El EH EH] does not occur. In the appendix [D] we employ the splitting (1.14) of the finite 



term to study the 0(1) term in (2.7) and discuss the approximation obtained by using the near 
horizon geometry. 



2.2 Warped black hole 

In this section we employ the observation just made about the role of the region close to the 



horizon and apply the expansion described in (2.1 ) and (2.2 ) to a black hole where only the near 
horizon geometry is known. 

In [25] a minimal consistent truncation of the type IIB supergravity has been considered by the 
following metric 

2 

dip + ^2 cos ' 



ds 2 



e-z x ds 2 M + R 2 e x 



-iii 
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\ 2 fi 2 



sin 



(2.10) 



where the non compact space M is given by 



ds 2 M 



ge 



"dt 2 + * 

9 



+ 



R 2 



(2.11) 



i=l 



The functions \i 9 an d w depend on the coordinate r only. The geometry (2.10) is required 
to provide AdS§ x T 1 ' 1 on the boundary, i.e. at large r. 

In [25] the equations of motion coming from the effective Lagrangian have been solved numer- 
ically; nevertheless analytic formulae have been found in some limits. We are interested in the 
T = regime, for which the first term of a series expansion near the horizon is given. The novel 
feature is that the near horizon region is a warped product AdSz x R 3 x T • . As discussed in 
|25j , one can employ the symmetries of the problem to set to one both the AdS radius and the 
position of the horizon, but we prefer to keep ro generic for clearness. 

The metric (2.10) falls into the general class considered in the section [l~T] through the ansatz 

9, d = 3 and 

A{r) 2 = g(r) B(r) 2 = * . (2.12) 



(1.16) by choosing D 



The analytic behavior near the horizon in the T = case reads [25 



9(r) 



b(r - r ) 13/3 + . . . 



b = 



93312 V12 
2J3 



(2.13) 



As checked in the section 2.1 for the charged black hole, the near horizon region determines 
the leading divergence of the finite term of the holographic entanglement entropy as the strip 
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in the boundary becomes large. Thus, we perform the expansion discussed at the beginning 
of the section [2] by using the near horizon geometry (2.13) instead of the full metric (which is 



still analytically unknown). Introducing r m \ n = tq + ep m \ a with finite p m \ n and changing the 
integration variable accordingly (r = r$ + ep), we get for the leading behavior of the integral in 



(1.20) the following result 



L 1 r° dp 1 0Fr(|) 

+ --- = ^7^wm^ ^ + ••■ ( 2 - 14 ) 



2 y/mVoe^ J Pmln p^(p - Pmin ) ' ' " VObTo^if ) (r min - r Q )^ 



where . . . denote higher orders in e. The same procedure can be applied to the integral in (1.23) 
which provides the leading divergence of the finite term in the holographic entanglement entropy 
as r m i n approaches the horizon. The result reads 



r 



5 \ 5/2 



-r(|) 



r\dr = -^= v yiJ + ... = rgL+ ... (2.15) 



yW)^^ -r^ iQ J V6bT(f ) (r min - r ) 5 / 3 



where in the last step we have used (2.14). Thus, also in this case the expected behavior for 



r m in — > tq is recovered (here we have d = 3). 



2.3 Perturbed Lifshitz background 

The Lifshitz background is defined by a metric which is scale invariant if the space coordinates 
and the time coordinate scale with a different power. The relative scale dimension of time and 
space is the dynamical exponent. This parameter usually occurs in the time component of the 
metric; therefore it does not affect the computation of the holographic entanglement entropy, 
which involves the metric on a constant time slice. An example of this type is considered in 
the section [3} Instead, when the dynamical exponent occurs in some spatial component of the 
metric, then it usually turns out to be involved non trivially in the holographic computation of 
the entanglement entropy |30j . In this section we consider an example of this type. 
A perturbation of the Lifshitz background through a formal parameter expansion was studied 
in [26] as a solution of the Abelian Higgs model in AdS^ [27] , introduced to describe supercon- 
ducting black holes. The metric to consider reads 

ds z = _ g{r) 2 dt 2 + f^ dS 2 + e 26(r) ^ ^2 ^ 
R z r z 

In [26] it was found that the Lifshitz background is a solution and also its perturbation of the 
following form is allowed 

9(r) = + A</i(r) + 0(A 2 ) b(r) = Xcr^ + 0(A 2 ) (2.17) 

where A is a formal expansion parameter and 7 depends on the dynamical exponent uj besides 
other parameters of the model. The explicit expression of 7 is not important for our discussion. 
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Notice that the dynamical exponent affects the spatial part of the metric through the perturba- 
tion of the Lifshitz background, and therefore it occurs in the computation of the holographic 



entanglement entropy. Since (2.16) on a constant time slice is a special case of the ansatz con- 



sidered in the section [lTJ we can employ the results discussed there. From (1.23) with d = 2, 
B(r) = r/R and A(r) = e~ b ^r/R, we find that the finite term in the holographic entanglement 
entropy is provided by the following integral 



,b(r) 



1 dr 



a/1 " (WO 4 



1 dr 



(2.18) 



+ Ac 



nn \A ~ ( r mm/r)- 



. dr + 0(A 2 



We are mainly interested in the 0(A) term in the r.h.s. of (2.18) because the O(l) one provides 
the result of AdS^ and of the Lifshitz background in four dimensions (they have the same 



entanglement entropy because their metric differs only in the time component). In (2.18) we 
cannot go to 0(A 2 ) because it involves the 0(A 2 ) of b(r) in (2.17), which is not known; but the 



O(A) term is already interesting because it contains the dynamical exponent through 7. To get 
a finite result from the integral at 0(A) in (2.18) when a — > 00 we need 7 < —1. Then we have 

-7 ^7+1 



J r„ 



in a/1 - (r m m/r)' 



dr 



Bp* 



1 + 7 1 
4 '2 



^r(^) 
4r(ip) 



„7+l 



+ 0(a 



1+71 



(2.19) 



where we found it useful to employ the integration variable p = r m [ Q /r and the final result is 
expressed in terms of the incomplete beta function B z (a, b), which reduces to the beta function 
B(a,b) = T(a)T(b)/T(a + b) for z = 1 and it is related to the hypergeometric function for a 
general z as B z (p, q) = (z p /p) 2 -Pi (p, 1 — q;l + p; z). 

As for the length L of the interval in the boundary, it is related to r m ; n through the integral 



(1.20), which in this case can be expanded up to O(A), similarly to what we have done in (2.18) 



for the area of the minimal surface. The result is 

Mr) 2 R 2 



L 



2R 2 r 2 ■ 

z, - rt 'mm 



r A x/ 1 - (rmin/r)' 



dr 



+ A 



Cj R 2 v^r(^) 



1-7t 
mm 



- 5—7 ' 

- 4 - 



+ 0(A 2 ) . (2.20) 



Again, the first term in (2.20) provides the result for AdS4 (see (A. 3)). We can invert (2.20) 
perturbatively and find r m j n (L) up to 0(A 2 ) terms by using that 



L 



» in [l + Cl ^ in A + 0(A 2 )] 



do 



A + 0(A 2 ) 



In our case we find 

_ 2R 2 ^T(l) 

rmin " ni)L 



1 + A 




4r(^)r(f) 



+ 0(A 2 



(2.21) 



(2.22) 
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Plugging this result into (2.19) we find that the correction O(A) to the holographic entanglement 
entropy is proportional to 



Ac 



un a/1 - (r m m/r)' 



: dr 



Ac V^T( 



-1— 



4r(i=^) 




l + 0(a 1+7 ) + 0(A 2 



(2.23) 

Since we are assuming 1 + 7 < this term diverges like L~( 1+ ^\ The interesting feature is 
that the dynamical exponent occurs in a non trivial way in the scaling of the finite term of the 
holographic entanglement entropy in terms of the width L of the strip. This computation is not 
conclusive because it involves only the first term of a perturbative expansion, but we expect the 
occurrence of the dynamical exponent in such scaling also for the result computed with the full 
(non perturbative) expression of the metric. 



3 A Lifshitz black hole in four dimensions 

In this section we consider the Lifshitz black hole in four dimensions (d = 2) found in 
Because of the simple emblacking function characterizing this black hole, we can compute the 
holographic entanglement entropy analytically to all order in the UV cutoff. This allows us also 
to check the method employed in the section [2] to find the divergent term in the finite integral 
of the area as z max goes to the horizon zq. 

The Lifshitz black hole of [2B] is a solution e.g. of a model in four dimensions which includes, 
besides gravity, a massive U(l) gauge field and a strongly coupled scalar, namely a scalar without 
kinetic term. Its metric reads 

. o _ , . dt> dx dz _ . . z , „ . 

The boundary is at z = and the range of the holographic coordinate is (0, zq). The dynamical 
exponent is ui = 2 and the bulk curvature radius R has been set to one. Near the boundary 



the metric (3.1) asymptotes the Lifshitz spacetime in four dimensions with dynamical exponent 



equal to two. Near the horizon the emblacking function f(z) vanishes linearly and the metric 



on the constant time slice is (1.3) with the f(z) given in (3.1). 

We remark that, since the anisotropy u) does not occur in the metric on the constant t slice, we 
do not see the effects described in [30 j. In that case they have an anisotropy between two spatial 
directions; therefore the holographic entanglement entropy is sensible to the difference between 
them. 



As first step we study the leading order for z max — > zq of the finite term (1.13) by employing the 
expansion described in the section [2} The result is 

A 2 (z mSLX , 0) = ^— log(z " ^max) + O(l) . (3.2) 

V2z 

Like in all the cases considered in the section [2j we cannot say anything about the finite term 
with this method. 
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For the Lifshitz black hole (3.1) we can compute the integral in (1.10) analytically (we find it 



convenient to adopt z = w /z^ax as integration variable). The result reads 



-42(>max,a) 



2z 2 

""max 



W 2 V 'fUw) \Amax - W 



: d,W 



I 



a 



where 



T(x) = 2^1 + PE\ arcsin 



2(1 -px) 



1 + 



(3.3) 



(3.4) 



1 + P 



' 1-x 2 
x(l — fix) 



being (3 = (z maiX / z$) 2 and the function F(x\m) and E(x\m) the incomplete elliptic integrals of 
the first and of the second kind respectively. Notice that the upper extremum of integration in 



(3.3) gives a vanishing contribution. Expanding (3.3) for small UV cutoff a we find 

3 



A 2 {z n 



2 

a 



MP) Pa 



z,. 



12 H 



a 



(-e 



43 
40 



,6 



(3.5) 



11 47 445 
28 56 448 / z 



a 



+ 0(a 9 



with the function fo(/3) occurring in the finite term of this expansion given by 



MP) = 2^/l+~pE ^arcsin 



V2[3K 



1 + P 



1 + P J V 2 

where K(z) is the complete elliptic integral of the first kind. As z max — > zq we get 

MP) 



V2 



log(z - z max ) + 0(1) 



(3.6) 



(3.7) 



which confirms the result (3.2) found through the method described in the section [2] 



For this black hole we can compute also the integral (1.8) as done for the one in (|3.3|) . Again, 



the upper extremum of the definite integral gives a vanishing contribution. The result reads 

2/3 



L 



2z„ 



PVi + P 



F ( arcsin 



yi+p 

V2 



1 + P 



(1-/3)11 



1 + P' 



; arcsin 



yr+P 

V2 



where II ( x, <^>|?7i) is the incomplete elliptic integral of the third kind. When z n 



L 



V2 



log(z 



Combining this result with (3.7) we obtain 



A 2 (z max ),0) 



L 



+ 0(1). 



+ 0(1) 



1 + P 
(3 

zo we have 
(3.9) 



(3.10) 



as expected. Besides providing another check for the method discussed in the section[2j this is the 
first case of a black hole whose holographic entanglement entropy can be computed analytically. 
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4 Two disconnected strips 



In this section we consider the case of a region A in the boundary made by two parallel strips. In 
particular, following [23], we study the transition of the mutual information in AdSd+2 (section 



4.1) and in the charge black hole background (section |4.2[ ). 

Let us consider a spatial slice of the boundary theory with two parallel strips A\ and A2 whose 
widths are L\ and L 2 respectively and separated by a distance Lq. As recalled in the intro- 
duction, the natural quantity to study for two disconnected regions is the mutual information 
Ma = Sa ± + Sa 2 — SA1UA2 because it is UV finite. 

In order to find the minimal surface associated to the region A = A\ U A2 in the holographic 
computation, we have to consider two pairs of disjoint surfaces extended in the bulk whose 
boundary coincides with the boundary of the two strips. Together with the region A, the first 
pair of surfaces encloses a connected volume of the bulk, while the second one encloses two dis- 
connected volumes of the bulk. The strong subadditivity inequalities guarantee that the pair of 
intersecting surfaces in the bulk whose boundary coincides with dA is not minimal |22 | 123 1 [23]. 
The divergent term giving the area law is the same for both these two pairs of surfaces because 
they share the same boundary. Thus, in order to find the pair with minimal surface, we have to 
consider the finite term (in the UV cutoff) of the integrals giving the area of the pair of surfaces. 



We find it useful here to change slightly the notation for the finite part (1.13) of the holographic 



entanglement entropy by introducing A d (L) = A d (z max , 0) where z max = z max (L) is the inverse 



function of (1.8). Thus, we consider 



S d (L 1 ,L 2 ;L ) = min[A d (L 1 ) + A d (L2);A d (L ) + A d (L 1 + Lo + L 2 )] (4.1) 

disconnected volumes connected volume 

which occurs in the mutual information for the finite parts 

M d (L x M-M) = A^) + A d {L 2 ) - SdiL^L^Lv) . (4.2) 



We remark that in (4.2 ) we talk about mutual information with a slight abuse of notation because 



the mutual information is given by \A.2\ multiplied by a factor i^L* -1 / '(4G^ +2) ) coming from 



(1.2) and (1.10). We made this choice for clearness and we believe it will not mislead the reader. 



The mutual information (4.2) is zero when the minimal surface is given by the pair of surfaces 
enclosing the disconnected volumes and it is positive when the minimal surface corresponds 
to the pair of surfaces enclosing the is the connected volume. The transition of the mutual 



information (4.2) from zero to a positive value occurs when the two terms compared in (4.1) are 
equal, i.e. 

A d {L x ) + A d (L 2 ) = A d (L ) + A d (L\ + L + L 2 ) . (4.3) 

In the remaining part of this section we study this equation in the special case of equal strips, 
namely L\ = L 2 . First we consider AdS d + 2 , where some analytic result can be found, and then 
the charged black hole in AdS d+2 . 
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Figure 3: Holographic mutual information in AdS± with Li = L\. On the left we show 
M 2 (Li,Li;Lo) for L = 0.87 (red), L = 1.91 (blue) and L = 3.93 (black). On the right, 
in the parameter space (L\, Lq), we plot the position of the transition point at which the mutual 
information starts to be non zero. 



4.1 AdS, 



d+2 



For AdS d+2 the analysis is simple because we explicitly know that (see ( |A.9 ) and (A.10)) 
ML) -- 



ad 
L d-i 



1 



ad 



V^r(^r 



d-\ 



r 



.2d) 



(4.4) 



which holds for d ^ 2. Keeping the distance Lq between the two equal strips fixed, for small L\ 
the pair of surfaces enclosing the disconnected volumes is minimal and the mutual information 



(4.2) is zero. Increasing L±, at a certain point the pair of surfaces enclosing the connected 
volume becomes minimal and the mutual information (4.2) is therefore positive. For large L\ 



the mutual information goes asymptotically to a constant, as shown for d = 2 in the figure [3] 
(plot on the left). In order to find the asymptotic value of the mutual information, we observe 
from (4.4) that Ad(L) — > when L — > oo. This implies that 



lim M d (Li,Li;L ) = -A d (L ) 



(4.5) 



which provides the asymptotic value of the mutual information as a function of the distance 
between the strips. 

As for the transition point at which the mutual information starts to be positive, its defining 
relation (4.3) specified for AdS d+2 and L2 = L\ reads 



(Lq/L, + 2 



(Lo/L^- 1 



2(L /L ] 



\d-l 



d > 2 



(4.6) 



where (4.4) has been employed. For any fixed d ^ 2, we can easily observe through a graphical 



analysis that the equation (4.6) has only one positive root for Lq/L\. This root provides the 



angular coefficient of the straight line in the plane (Li,Lo). In the figure [3] (plot on the right) 
the case of AdS^ is considered. 
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Figure 4: Two equal and parallel strips. Angular coefficient of the line characterizing the transi- 
tion of the holographic mutual information in AdS d +2 i n terms of d. The red point corresponds 



to AdSs, which is not described by the equation (4.6). In this case the transition occurs at the 
value x = 1/2 of the conformal ratio 



In the figure [4] we show the angular coefficient of the straight line, namely the solution of (4.6), 



as function of d. We remark that the equation (4.6) holds for d 2. The case of AdSs (i.e. 
d = 1) has been studied in |24| . finding that the transition occurs when the conformal ratio 
x = zi2Z34/(zi3Z2i) = L\/{L\ + Lq) 2 = 1/2, which corresponds to the red point in the figure|4j 

4.2 Charged black holes 

In this section we consider the holographic mutual information for a charged black hole. 
By employing the results of the section [2j we have that 

A d (L) = ^+c d + o(l) for large L (4.7) 



where c d is the 0(1) term in (2.5), (2.7) and (2.9). Since we are not able to determine Cd 
analytically, we fix it by fitting the numerical values of A d (L) at large L with a line. 
For two equal strips of width L\ at fixed distance Lq, the behavior of the mutual information 
is qualitatively the same obtained for AdSd+2 and shown in the figure [3] (plot on the left). The 



asymptotic value of M d (Li, L±; Lq) at fixed Lq can be found by employing (4.7). For L\ — >• oo 
we have 

M d (L u Li; Lq) = 2A d (Lx) - A d (L ) - A d (2L 1 + Lq) — ► c d - A d (L ) - ^ . (4.8) 

As for the position of the transition point of M d (Lx, L\\ Lq) in the plane (L\,Lq), the curve is 
instead qualitatively different from the corresponding one obtained for AdS d+ 2, Indeed, while 
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Figure 5: Extremal charged black hole in AdS^ with zq = 1. On the left the position of the 
transition point of the holographic mutual information in the parameter space (Li, Lq). On the 



right, a zoom of the same plot: the asymptotic line is provided by the equation (4.9) and the 
green line corresponds to the transition point of AdS^ (figure [3j plot on the right). 

we get a straight line for AdSd+2 (plot on the right in the figure [3]) , for the charged black hole we 
find a curve with an asymptotic constant value (plot on the left in the figure [5]). In particular, 
the straight line of AdSd+2 is tangent to the curve corresponding to the charged black hole which 
is asymptotically AdSd+2-, as shown by the plot on the right in the figure [5j Indeed, for small 
values of L\ the pairs of surfaces to compare are close to the boundary and consequently the 
transition between them is determined by the asymptotic geometry. 

Let us consider further the characteristic asymptotic value Lq of the curve of the transition 
points of the mutual information for a charged black hole in the plane (L\,Lq) as L\ becomes 
large. The equation defining Lq can be found by taking the limit L\ — >• oo and Lq — > Lq of the 



equation (4.3) and employing (4.7). The result is 



MLo) + - d 



(4.9) 



which can be solved numerically. This asymptotic value of the distance between the two strips 
could be interpreted as a signal of the occurrence of a finite correlation length in the boundary 
theory. 

The qualitative features just described for the extremal charged black are found for the non 
extremal case as well. The mutual information Md(Li, L\\ Lq) behaves like in the plot on the 
left of the figure [3] and the curve of the transition points is qualitatively like the one shown 
in the figure [5j with the asymptotic value given by the solution of the equation (4.9) with the 
proper emblacking function depending on the temperature. In the figure [6] we show the curves 
of transition points of M^{L\, L\\ Lq) for two different temperatures besides the extremal case 
at fixed charge. The curve corresponding to a certain temperature always stays below the curve 



corresponding to a lower temperature, meaning that the asymptotic value determined by (4.9) 



decreases with the temperature for a fixed charge of the black hole. We recall that imposing Q 
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Figure 6: Charged black hole in AdS± at fixed charge. Position of the transition point for 
M2(.Li, L±; Lq) in the plane (Li,Lo): the red curve corresponds to the extremal case (see the 
figure [5j plot on the right) while the blue and the orange ones correspond to two non extremal 
cases (T = 0.1 and T = 0.18 respectively). The black horizontal line gives the asymptotic value 
of the extremal case while the green line corresponds to the transition point of AdS^. 



5 Conclusions 

In this paper we have considered two aspects of the holographic entanglement entropy in black 
hole backgrounds: the behavior of the finite term as the width of the strip in the boundary 
becomes large and the transition of the mutual information for two equal strips in the parameter 
space given by the width of the strips and their distance. 

For one strip in the limit of large volume, which means that the turning point of the minimal 
surface approaches the horizon, we confirm and extend to new cases the known result that the 
finite term scales like the width, and therefore like the volume, of the strip. The distinguished 
feature of the different black holes is the degree of the divergence of the finite part in terms of 
the distance between the turning point and the horizon, which is determined by the near horizon 
geometry. In the case of a Lifshitz background with a dynamical exponent entering in the spatial 
part of the metric, such scaling could be influenced by this exponent. For a Lifshitz black hole in 
four dimensions we computed the analytic expression of the holographic entanglement entropy 
to all orders in the UV cutoff. 

For two equal and parallel strips in the boundary, we have found that the transition of the mutual 
information for a charged black hole naturally provides a finite limiting distance between the 
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strips as their width becomes large. This asymptotic value could be interpreted as a signal of 
a finite correlation length in the boundary theory. The transition in the mutual information is 
characteristic of the holographic prescription; therefore it is a large c effect. We believe that 
it is important to further study this transition in order to understand how it smooths out for 
finite c. This is part of the general aim of reproducing through holography the results obtained 
for the mutual information in the finite c models. 

Acknowledgments 

It is a pleasure to thank Hong Liu for collaboration in the initial stage of this project and for 
many helpful discussions and advices during its development. I am also grateful to Matthew 
Headrick, Mark Hertzberg, Veronika Hubeny, John McGreevy, Mukund Rangamani, Tadashi 
Takayanagi, Frank Wilczek and in particular to Antonello Scardicchio for useful discussions. I 
thank the Physics department of the University of Pisa and the Galileo Galilei Institute for their 
kind hospitality during the last part of this project. 

This work is mainly supported by Istituto Nazionale di Fisica Nucleare (INFN) through a Bruno 
Rossi fellowship and also by funds of the U.S. Department of Energy (DoE) under the cooperative 
research agreement DE-FG02-05ER41360. 



A AdS d+2 

For the sake of completeness, in this appendix we briefly review the results for the holographic 
entanglement entropy in AdS^+2 f° r the strip [BJ [7]. The expressions in the section [I] can be 
applied with f(z) = 1 identically. 

The inverse function of the profile z(x) representing the minimal surface is given by 



, w d , w d+1 „ fd + l 1 3d + l w 2d 

x ( z ) = / -^F=^ri dw = 71 , Tv^ 2^1 



y/z^-w™ (d+l)zi ax * x \ 2d '2' 2d 'zZ x 

r(^) z d+1 fd + i l sd + i z 



(A.l) 



r(&) (d + l)zL x 2Fl { 2d '2 ; 2d ; zggj ■ ,A/2) 



Since x(0) = L/2, from (A. 2) we see that 



L — ( J_\ d ^ 2max • (A-3) 

1 \2d) 



As for the regularized area of this minimal surface, it is given by (1.10) where now the integral 
to perform is 



•Smax ~d 



Mz ma „a) = 2/ ™„dw (A.4) 



V 



a W^x/Z^-W 



+ i a I ,„ - 1 I dw (A.5) 



(d - 1) a*"* (d - 1) ztJc J a W d \ _ W 2d 
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where the divergence for small a has been isolated as in (1.11) or ( 1.14[ ) (in absence of the black 
hole they provide the same result). The integral in (A. 5) reads 



f 

■J a 



V 



*max w 



1 dw 



(A.6) 



(d- ljto 4 *- 1 
2 



1 - 4/1 



2d 



2w 



d+l 



y 2d 



id? l)^ma: 



2^1 



d+l 1 _ 3d + 1 10 



2f/ 



2d '2' 2d 's£L 



-j — r ™' . = + O a 

d-1^ d-i)r(i ^ 



<M-1n 



d ^ 2 . 



(A.7) 



Thus the UV divergence of Area^^) has been isolated and the final result is [7] 

2 



(d- i) a*-' (<i-i)r(i)4ix 



2 ^ r <^) r + o(</«) 



(d - 1) c^- 1 L d - 



T + 0(a' 



where 



r(i) 



d- 1 



(A.8) 
(A.9) 

(A.10) 



This expression has been employed in the section [471] to study the asymptotic value of the mutual 
information. 



Now we find it useful to derive (A.8) also in the following way, which could be employed in a 



generalized version for the black holes. First one writes the integral in (A.5) as a series 

lax 1 / \ 



1 dw 



E2b n 
7 2dn 
n=l max J a 



w^^dw 



'max 



(2n - l)d+ 1 

n=l 



+ 0(a 



(A.ll) 
(A.12) 



where in ( A.ll ) the coefficients b n can be found by employing the following identity with a = 1/2 



1 



(i- x y 



2 F 1 (a,/3;(3;x) = 



(a)n 



X 



(l/2)n 



n=0 



111 



111 



(A.13) 



being (c) n = c(c + 1) . . . (c + n — 1) the Pochhammer symbol (we recall that (c)o = 1). 
Then, from (A.5) and (A.12) we get that for the finite term in the expansion for small a 

2 / 



y d-l \ A _ i + ^2 



^ (2n-l)d + l 

n=l , 



(d-i)r(i)^ 



(A.14) 



which agrees with the finite term in (A.8). 
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B Charged black holes in AdSd+2 



In this appendix we review some features of the charged black holes which are asymptotically 
AdSd+2- The metric reads 

,2 r 2 , , IHaN R 2 dr 2 , Q 2 M ,„ „, 

^ = ^(-M 2 + ^ 2 ) + ^ T / = 1 + ^- r ^T (B.1) 

for d ^ 2, where dx 2 is the metric of M. d , M is the mass and Q is the charge of the black hole. 
The boundary corresponds to large r, where the metric becomes the one of AdSd+2 with radius 
R. The Schwarzschild black hole in AdSd+2 is obtained by setting Q = 0. 

By introducing the variable z = R 2 /r, the metric (B.l) becomes (2.3) and the boundary corre- 
sponds to z = 0. This parameterization is largely used in this paper. Another useful parame- 
terization of the radial coordinate is 

e" /R = £ ■ (B.2) 
Notice that a scaling of r corresponds to a shift of f. With the parameterization given by f, the 



metric (B.l ) reads 

ds 2 = e f ' R [-mdt 2 + d^] + ^ f(r) = l + ^e- 2 ^__^_ e -(. + i)^. (b.3) 
It is convenient to parameterize Q by introducing follows 



„ d+1 2d d + 1 (R< 



2\ 2d 



From this expression it is evident that Q has the dimension of [L] d . The limit z* — > oo corre- 
sponds to the Schwarzschild black hole. The chemical potential reads 



^ = / d g F Q _ v / d(rf + 1) g F r Q f r *\ d _ V d (d + 1) qf f z °\ d (35) 



2{d-l) i?2 r d-i ■ v^(rf-l) R2 VoJ V2(d-l) z 



'o 

where g F is the effective dimensionless gauge coupling. When z* — > 00 for fixed zo the chemical 
potential vanishes. The temperature is 

T= (d + l)r ( _r?\ =d± l( d+l/ _d-lQ^\ 

AttR 2 V r 2 V 4ttz V 4 d ) 4vrz V d+1 ^ [ ' 

Since ro r* in order to impose T > 0, we have zq ^ z*. Notice that if we want to keep Q 
fixed, changing T implies a change of Zq. Indeed, the values of Q and T fix the position zq of 



the horizon through (B.6), which can be written also as follows 

(d-l)Q 2 2d 4vrT 
(dTTW^ d+1 20 ~ 1 ~ ' 

Setting R = 1, if we decide to choose zo = 1 at T = then Q 2 = (d + l)/(d — 1). Keeping this 
value for Q 2 fixed, moving to T > modifies zq according to (B.6) which becomes 

z o" + !f^o-l = 0. (B.8) 
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From the relation (B.6) it seems that there is a maximum temperature corresponding to —¥ oo. 

V2(d+l)(d-l) 



Instead the relevant parameter is the ratio 



4 d \ 

?2d I 



(B.9) 

which spans all the positive real numbers when £ [jzq, oo) in a strictly monotonical way, going 
to infinity when — > oo. From (B.9) we can see that (the other root is negative) 

d 



T 



2a d [i 



+ 



T 



2a d fi 



+ 1 > 1 



(B.10) 



which becomes 1 when T = for any d ^ 2. The parameter M, which can be expressed in terms 
of Q 2 and the position tq of the horizon, reads 

2d / r>2 \ d +! 



M 







+ 



„d-l 



r d+l 

r o 
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d+1 r. 



d-lrt 1 



'0 " " '0 

Thus, the emblacking function can be written as follows 
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(B.12) 
(B.13) 



Notice that from (B.12) and (B.10 ) we can write f(z) in terms of the ratio T/ (jl. A very important 
role in our discussions is recovered by the near horizon geometry, namely the one obtained when 
z — > zq. Close to the horizon, the emblacking function can be expanded as 



/(*) 



(d + DW - 4% _ Zo) + ^ + _ 2o)2 + 0((2 _ Zo)3) (B 14) 



Zq 



2z 2 



4irT(z -z) + 



d(d + 1 - 6ttzqT) 



(z-z ) 2 + O((z-z ) 3 ) • 



(B.15) 



In the extremal case (T = 44> z* = zo) the emblacking function /(z) = 0((z — zq) 2 ^, while in 
the non extremal case (T > and z* > zq) we have /(z) = 0{z — zq) when z — > zq. We remark 
that also in the case of the Schwarzschild black hole, which corresponds to z» — > oo, we have 

f(z) = 0(Z - Zq) as Z-> Zq. 



C Disk geometry 

In this appendix we briefly discuss the case in which the region A in the spatial section of the 
boundary theory is given by a disk, while in the bulk a black hole occurs whose metric on the 



constant time slice is given by (1.3). 
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Taking as A the circle given by p = R (it is more convenient to adopt the polar coordinates) 
and assuming that z = z(p), we get 



Area(7^ 



V d . 



dpp 



d-l 




(C.l) 



where z' = dz/dp and Vd-\ is the volume of the d — 1 unit sphere. Now the Lagrangian density 
£disk[<z(p)] is the integrand of (C.l) and it explicitly depends on the coordinate p. This means 
that there is not a conserved first integral. 



In order to minimize the functional (C.l) we need to solve the second order equation given by 
the equation of motion, which is 
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+ {d-i) 



z z 



1 + 
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+ dp 



1 + 



p vW f{z) 







(C.2) 



where f'{z) = df(z)/dz. Thus, this case is more complicated than the strip, largely considered 
throughout the paper, because now we have to solve a second order equation to find the profile 
to use in the integral giving the area. 

For AdSd+2 the equation to solve is (C.2) with f(z) = 1 identically (see the footnote 20 of |7J) 
and its solution reads 

z (p) = V^-P 2 (C3) 
which is the semispherical surface whose A is the maximal circle. For a black hole background, 



which has a non trivial emblacking function f{z) , the equation ( C.2 ) for the profile of the minimal 
surface can be solved numerically. 



D An alternative splitting of the finite term 



In this appendix we provide some insights about the expansion for z max — > zq of the finite term 
of the holographic entanglement entropy and about the role of the near horizon geometry by 



considering the splitting (1.14). 



Let us assume to know the first integral in ( |1.15[ ) analytically. Then, the O(l) term of ^4f/(^maxi ^0 
in the expansion for a — > is obtained by ld(0, z max ) plus a contribution from the first integral. 
In general we are unable to compute ld(0, z max ). Anyway, we are interested into its expansion 
as z max — > zq. The emblacking function f(w) depends on the ratio w/zq. By introducing 
y = w/z max £ [0, 1] as integration variable, the function /(z max y) depends on the ratio z max /zQ < 
1, therefore we can consider the expansion of the function 1/y f(z max y) as z max / zq 
obtaining 
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(D.2) 
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Unfortunately, the integral and the series cannot be inverted because the integrals occurring for 
any fixed n are divergent at the upper extremum y = 1 as we will see below in a special case. 
By introducing an intermediate scale a < z\ < 3 max , we can write 



(D.3) 



Now, in 1(0, z\) we can invert the series and the integral because the upper limit is z\/z ina _ x < 1 
and the integrals converge. We get 
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1 dy 
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30 



(D.4) 



which is a well defined expansion whose coefficients depend on the ratio 3^/z max . 

The second integral I(z\,z max ) is still divergent when 3! max /3o — > 1 _ and we cannot invert the 

series with the integration as done in (D.4); therefore it must be computed analytically. Since 

this is usually too difficult, we can approximate it by employing the near horizon behavior of 

the emblacking function. The closer is z\ to 3 max , the better is this approximation. 

In order to apply these considerations to a concrete example, let us consider the extremal charged 



black hole in AdS^. The first integral in (1.15) in this case can be computed, obtaining 

"^max Q 
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(D.5) 



log 



Aw + 2z + y/6(3w 2 + 2z w + zj) \ ^3w 2 + 2z w + 



+ 



\/6z 



log(3 



ZQ-W 



+ 0(1) 



30 W 



+ 0(a 2 



(D.6) 



where the square brackets in (D.6) enclose the finite term in the power series in a, which has 
been further expanded for z max — > zq. 

Now, by expanding the integral in (1.15) as explained in the section [2] we find 

7T 2 

h(Zmax, 0) = 7 = . + log(z - 3; m ax) + 0(1) 

VO3 y/Z Q - Z mSuX V6^0 



(D.7) 



where, again, we do not control the finite term. Notice that the logarithmic divergence in (D.7) 



cancels the one in (D.6) and the remaining divergence is the same one found in (2.7) by using 
( |1.11[ ). This is a consistency check of the two splittings ( |1.11[ ) and (1.14) of the same integral. 
As discussed above in this appendix, let us consider the integral 12(0, z m3iX ) in terms of the 
variable y (see (D.l)). The emblacking function then reads 

3 



/(^max U) 
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30 
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30 



(D. 



By expanding 1/a/ /(3 max y) for 3 max / zq — >• 1 , we find the functions h n (y) occurring in the 



series (D.4). For the first terms, they are e.g. 
ho(y) 



v/1 - 4y 3 + 3y 4 



hi(y) 



6(1 -y)y 3 



(1 - 4y 3 + 3y 4 ) 3 / 2 



(D.9) 
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and the corresponding integrals obtained by inverting the summation and the integration in 
(D.2) are divergent in 1 because 1 — 4y 3 + 3y 4 = 0((1 — y) 2 ) when y 



1. 



As discussed above, we introduce an intermediate scale Z\ and split the integral as in (D.3), 
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Figure 7: Extremal charged black hole in AdS± and zq = 1. Plot of l/v/ f(z) (solid line) and 
W /nhC-z) (dashed line) as functions of z G (0, zq). 



obtaining for the first term a well defined power series (D.4) in terms of integrals involving the 
functions h n . We are not able to compute them analytically, but we are guaranteed that in 
-^2(0, z\) is finite as z max —> zq. The divergence comes from the near horizon region. 



As for the second integral in (D.3) giving the divergent part for z max — > zq, we cannot compute it 
explicitly, but we can relate it to the corresponding integral involving the near horizon geometry. 
In particular, as shown in the figure [7J the integral hizx, z max ) is greater than the corresponding 
one computed with the near horizon geometry for any choice of zy, namely 

rzm&x 2 1 z 2 \ 

h(zx,Z mSkX ) > / 2) nh(^, *max) = / ~ r-z — -. — F , . = . ~ 1 ) dw (D.10) 

Jzy W Z \/_tr,h(w) \ \/zi„ - W 4 / 



where the emblacking function close to the horizon reads (see (B.15)) 



6 (w z ) 2 



Z 



The integral in (D.10) is easier to deal with and the closer z\ is to z max the better is the 



approximation obtained by substituting hizx, z m3bX ) with /2,nh(^Aj -Zmax)- 
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